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Abstract 

The blackbody is considered in the external general field. The additional coefficients of 
stimulated emission and stimulated absorption are introduced into the Einstein mechanism 
and the generalized Planck formula is derived. The Einstein and Debye formula for the 
specific heat is possible to generalize. The Bose-Einstein statistics is broken in the external 
field. The relation of the theory to the sonoluminescence, the relic radiation and solar 
spectrum is considered. 

Key words. Planck formula, Einstein's blackbody, phonons, Einstein's and Debye's 
specific heat, sonoluminescence, relic radiation. 



1 Introduction 

The distribution law of photons inside of the so called blackbody was derived in 1900 by 
Planck (1900, 1901), (Schopf, 1978). The derivation was based on the investigation of the 
statistics of the system of oscillators. Later Einstein (1917) derived the Planck formula 
from the Bohr model of atom. Bohr created two postulates which define the model of 
atom: 1. every atom can exist in the discrete series of states in which electrons do not 
radiate even if they are moving at acceleration (the postulate of the stationary states), 2. 
transiting electron from the one stationary state to other, emits the energy according to 
the law 
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TlL) = E m - E n , (1) 

where E m is the energy of an electron in the initial state, and E n is the energy of the final 
state of an electron to which the transition is made and E m > E n . 

Einstein introduced coefficients of spontaneous and stimulated emission A mn , B mn , B nm . 
In the case of spontaneous emission, the excited atomic state decays without external stim- 
ulus as an analogue of the natural radioactivity decay. The energy of the emitted photon 
is given by the Bohr formula (1). In the process of the stimulated emission the atom is 
induced by the external stimulus to make the same transition. The external stimulus is a 
blackbody photon that has an energy given by the Bohr formula (1). 

If the number of the excited atoms is equal to N m , the emission energy per unit time 
conditioned by the spontaneous transition from energy level E m to energy level E m is 

Pspont. emiss. N rn A. rnr fiuJ , (2) 

where A mn is the coefficient of the spontaneous emission. 

In case of the stimulated emission, the coefficient B mn corresponds to the transition of 
an electron from energy level E m to energy level E n and coefficient B nm corresponds to 
the transition of an electron from energy level E n to energy level E m . So, for the energy 
of the stimulated emission per unit time we have two formulas : 

Pstimul. emiss. Qu)N m B mn h,U) (3) 
Pstimul. absorption Qui N n B nm huJ . (4) 

If the blackbody is in thermal equilibrium, then the number of transitions from E m to 
E n is the same as from E n to E m and we write: 

N rn A mn hu + N m Q w B mn hu) = N n Q u B nm tiW, (5) 

where is the density of the photon energy of the blackbody. 
Then, using the Maxwell statistics 

N n = De~&, N m = De-^, (6) 

we get: 



Qu, = B B Z • (7) 

The spectral distribution of the blackbody does not depend on the specific atomic 
composition of the blackbody and it means the formula (7) must be so called the Planck 
formula: 

hm 3 1_ 

vl C gfeT 



After comparison of eq. (7) with eq. (8) we get: 
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7T 2 C 3 

B m n B nm foijj-i ^~ rnn ' (^) 

It means that the probabilities of the stimulated transitions from E m to E n and from 
E n to E m are proportional to the probability of the spontaneous transition A mn . So, it 
is sufficient to determine only one of the coefficient in the description of the radiation of 
atoms. 

The Planck law (8) can be also written as 

q u = G{u) < E u >= U ^ . (10) 

ix 2 c A ewr - 1 

where the term 

<E U >=^- (11) 
ekr — l 

is the average energy of photons in the blackbody and 

G(u) = 4-s (12) 

is the number of electromagnetic modes in the interval u, u + du. 

Let us remark that coefficients A mn of the so called spontaneous emission cannot be 
specified in the framework of the classical thermodynamics, or, statistical physics. They 
can be determined only by the methods of quantum electrodynamics as the consequences 
of the so called radiative corrections. So, the radiative corrections are hidden external 
stimulus, which explains the spontaneous emission. 



2 N-dimensional blackbody 

The problem of the N-dimensional blackbody is related to the dimensionality of space, 
or, of space-time and some ideas on the dimensionality was also expressed and analyzed 
in recent time by author (Pardy, 2005). 

The experimental facts following from QED experiments, galaxy formation and 
formation of the molecules DNA, prove that the external space is 3-dimensional. With 
regard to the Russell philosophy of mathematics, there is no possibility to prove the 
dimensionality of space, or, space-time by means of the pure mathematics, because the 
statements of mathematics are nonexistential. The existence of the external world cannot 
be also proved by pure mathematics. However, if there is an axiomatical system related 
to the external world and reflecting correctly the external world, then, it is possible to 
do many predictions on the external world by pure logic. This is the substance of exact 
sciences. We know for instance that the success of special theory of relativity is is based 
on the adequate axiomatical system and on logic. 

In case of the n-dimensional blackbody, the number of modes can be determined (Al- 
Jaber, 2003). We use here alternative and elementary derivation. We know, that in case 
that the electromagnetic field is in a box of the volume L n , the wave vector k is quantized 
and the elementary volume in the k-space is A 0n = (2n) n /L n . 

The elementary volume of the n-dimensional k-space is evidently the volume dV n 
between spheres with radius k and k + dk: 
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2W 2 1 2vr n / 2 
dV n = d { ^k n } = —r-^k^dk, (13) 

Mi) J r (t) 

where T(n) is so called Euler gamma-function defined in the internet mathematics 
(http: / / mathworld.wolfram.com/GammaFunction.html) as 
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r(x) = J q e-H^dt- r(n/2) = (W 2(w y . (14) 

The number of electromagnetic modes involved inside the spheres between k and k + dk 
is then, with uj = ck, or k = uj/c and dk = dcu/c, 

2x^ = 2x ' ' ;2 L-^, (15) 

A „ r (f ) 2( n - 1 )W 2 c n y J 

where isolated number 2 expresses the fact that light has 2 polarizations. It means that 
the number of modes G n , in volume V, involved between spheres k and k + dk, or, uj and 
uj + cL;, is 

111 cj n_1 
G n (o;)dc> = 2 x p — —r — j — r — -L n duj. (16) 

For the spectral distribution, we have 

, _. G(uj)duj 

g nuJ duj =< E n > — — — =< E n > g{uJ)duj, (17) 

where the average value of energy < E n > is given by eq. (11). 

The spectral distribution of the n-dimensional blackbody has the final form following 
from eq. (17): 

hu I 1 1 uj 11 " 1 
Q ™ ~ 2(- 2 ) r (f)W 2 c- ' ( } 

3 Blackbody in the external field 

The Einstein derivation does not consider the situation where the blackbody is influenced 
by some external field. Einstein was not motivated by the case of so called sonolumi- 
nescence, or by relic radiation and so on, which play at present time in physics the 
fundamental role. Now, if we include some external nonspecified field, it is necessary to 
introduce transition coefficients C mn , C nm ,. From the equilibrium condition 

■N m -A mn -\- N m Q L0 B mn -\- N m C mn N n Q UJ B nm -\- N n C nm (19) 
and Maxwell statistics (6) we get: 

A-mn + QujB mn + C mn = e kT Qu)B nm + e kT C nm , (20) 
which can be modified using the definitions 
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ir 2 r 3 r r 

r> _ r> _ !__ a ■ mn — a ■ nm —ft (o\\ 

J - > nm J-'mn ft 3 mn i a Hmm * Hnm v^ 1 / 

as follows: 

1+Pmn- e^finm = 6"J^ \ ekT ~ l ) ■ ( 22 ) 

The generalized Planck law follows from the equation (22): 

^ " ^ e |f _ 1 + ^3 e W - 1 ' 1 J 

where P = j3 mn , Q = j3 nm are some function which must be calculated using the advanced 
solid state physics and advanced quantum mechanics. In order to get the second term in 
eq. (23) finite for u — > oo, it is possible to postulate the mathematical structure of Q(ou) 
as follows: Q{uj) = Rexp(-Su), where R, S are some constants. 

4 Discussion 

The original Planck derivation of the blackbody radiation was based on the relation 
between the entropy of the system and the internal energy of the blackbody denoted by 
Planck as U. 

While from the postulation of the relation 

d 2 S const , jN 

MP=-— (24 > 

the Wien law follows, the a priory generalization of eq. (24) gives new physics. The 

generalization was postulated by Planck in the following form: 



d 2 S k 



dU 2 U(e + U) ' 
The first integration can be performed using the integral 



I x(a 



dx 1 , 
= In 



;(a + bx) a 

We get the following result 



- + b 

x 



(25) 



(26) 



T dU e \U 

The solution of eq. (27) is 



1 dS k , / e 

In - + 1 . (27) 



e s/kT _ l 

The general validity of the Wien law 



U = -7^^ (28) 



£ - -/ f" ) (^) 



dU v \v 
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confronted with the equation (27) gives the famous Planck formula 



e = hu. (30) 

The next step of Planck was to find the appropriate physical statistical system which 
led to the correct power spectrum of the blackbody. This model was the thermal reservoir 
of the independent electromagnetic oscillators with the discrete energies. 

We know that Einstein later identified the reservoir with the oscillators in the 
crystalline medium in order to get his famous heat capacity with the high temperature 
limit via the Dulong-Petit law. He supposed that the frequency of the crystalline medium 
is ujei &e being the Einstein frequency. Then U = 3N < E >, N being the number of 
oscillators. Then, 

dU 3Nk(^) 2 eM^ E /kT) 
V dT (exp(hu E /KT) - l) 2 ' 
For kT >> Tiuje it is exp(hujE/kT) m 1 + fiuiE/kT and we get 

C v « 3Nk. (32) 
For kT << Tiuje we get from the formula U the following formula the approximation 



U « 3Nhu E e~-^, (33) 
from which 

dU 3Nkh 2 uj 2 E 

Cv = W = ^rf e (34) 

and evidently Cy — > 0, when T — > 0, which is in a agreement with experiment. 

However, the absolute harmony of theory with experiment was not achieved, because 
experimentally it is C v ~ T 3 , which does not follow from the formula (34). 

Later Debye in 1912 generalized the Einstein model in order to get more realistic 
formula for the heat capacity of every crystalline matter (Reissland, 1973). Debye 
assumed that the solid state can be represented by the continual medium with many 
modes. However, if we consider a continual medium, then there are the infinite number 
of frequencies. Debye resolved this contradiction by the normalization condition that the 
number of frequencies must be 3N. Or, 

/ G(uj)duj = 37V, (35) 
Jo 

where G{uj) is the number of phonon modes in the three-dimensional continuum. Or, 

x , V ( 2 1 \ 2 , 3tu 2 V , 
G(u))du) = — - —5- H — 5- duj — — (36) 
V ; 2tt 2 \vf vf J 2tt 2 c 3 ' V ; 

where Vt is the transversal velocity with two polarization and vi is the longitudinal velocity 
of sound and we have defined c by the equation 
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? = U + ?J- (37) 



From equation (35) follows the Debye maximal frequency 



9 N ( 2 1 . 



The average energy is then of the form: 



U = 



r^r, (39) 
Jo e~ — 1 



2tt 2 c 3 

which can be written using the substitution x = huj/kT in the following form: 

J^WpA (40) 
2wVh 3 Jo e'-V [ > 

The natural approximation is for kT << hu^ where the upper limit of the integral is 
infinity and we have with 



00 X 3 dx _ 7T 4 

o e x - 1 ~ 15 ^ ' 

u = aTi -- a = ^ff v - (42 > 

The corresponding specific heat in the last approximation is Cy = dU/dT = AaT 3 . 
It can be easily to show that if we introduce the debye temperature by the relation 

T D = ^, 0* = ^— J (43) 

then 

The corresponding heat capacity is then given by the formula 

C v = 3Nk{ D (^)-^D>(^)}. (45) 

The Debye formula is only an approximation because the frequency spectrum of the 
crystalline and noncrystalline medium is very complex. There are also many defects in 
the crystalline matter and it means that the determination of this spectrum is not easy 
mathematical problem. Some theory of the determination of such spectrum was discussed 
for instance by Maradudin (1966) and others. 

Nevertheless, after some elementary approximation, Debye was able to derive the heat 
capacity of his model of crystal in the following form (Rohlf, 1994): 

C v = ^Nk T D = hu D /k. (46) 



D 
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Of course, the Debye model is in the better agreement with experiment that the 
Einstein model. 

The formula (23) can be easily written in the form 



Qu, = 2 3 - (1 + P(w) - Q(w)ei*), (47) 

7T Z C° e ~kT — 1 

which means that the spectrum of the electromagnetical modes inside of the blackbody is 
modified. If we want to involve the external field in the Einstein model, Debye model and 
more sophisticated models, then the formal operation consists in the following elementary 
transformation 



CkT — 1 ekT 



-(l + P(u)-Q(u)e&) (48) 



If this formal procedure has the physical meaning, then it is the generalization of the 
Bose-Einstein statistics (BE). The known generalization of the BE statistics and Fermi- 
Dirac (FD) statistics is known as the parastatistics. It is defined as the statistics where 
in one energetical state can be p particles. If p is arbitrary, we get the BE statistics and 
if p — 1, we get so called FD statistics. The derivation of parastatistics is involved in 
physics as an exercise of combinatorics (Isihara, 1971). 

Here, we consider the statistics which is generated by dynamics. The historical 
analogue of this approach is the Maxwell statistics based on the dynamics of particles 
in reservoir and Einstein statistics (BE) of photons inside of the blackbody, derived on 
the basis of the dynamics of emission and absorption of photons. Our derivation is also 
based on the dynamics of emission and absorption of photons with the additional terms. 

The different generalization of the Bose-Einstein statistics was realized for instance by 
Bekenstein (Bekenstein et al., 1994) in order to explain so called the gray-body radiation. 

However, the idea of modification of the Bose-Einstein statistics or is forbidden by the 
quantum field theory where there is the strong connection between spin and statistics. 
The particles with the integer spins are bosons and the particles with the half-integer spins 
are fermions (Berestetzkii et al., 1999). Only if we accept the idea that the statistics of 
photons and phonons is of the dynamical origin, then we can also write the generalized 
heat capacity of solid state. 

The next possibility to modify the statistics is to write the Planck hypothesis (25) in 
the generalized form 

d 2 S _ a 

dlP = E n a n U n { } 

and then to determine the modified Planck law for the specific thermal systems. To our 
knowledge this way was not still realized. 

The formula (23) can be also related to the solar spectrum which is according to the 
NASA expertise not absolutely identical with the spectrum of the blackbody (Thekaekara, 
1994). 

The formulas for the modified specific heat can be eventually considered in case of 
the investigation of the specific heat of the carbon nanostructures in some nonspecified 
external field. Some investigation was performed by Li and Chou (2005). The theory 
and experiment with the carbon nanonstructures is at present time very actual area of 
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investigation and it means that every new information of the physical properties of such 
structures is useful. 

The formula (23) can be probably related to the spectrum of the sonoluminescence, 
because the external field is the ultrasound. Sonoluminescence is defined as the emission 
of short burst of light from imploding bubbles in a liquid when excited by ultrasound. The 
effect of sonoluminescence was discovered at the University of Cologne in 1934. Franzel 
and Schultes put an ultrasound transducer in a tank of photographic developer fluid. 
They hoped to speed up the development process. Instead, they noticed tiny dots on the 
film after developing, and realized that the bubbles in the fluid were emitting light with 
the ultrasound turned on. The bubbles are very small when they emit the light. About 1 
micrometer in diameter. The high compression of a small bubble of fluid is similar to the 
explosive compression of a pellet of material by laser beams, one of the methods proposed 
for the nuclear fusion. 

The situation in cosmology is a such the "cosmic microwave background" (CMB) 
discovered in 1965 by Penzias and Wilson is of the Planck distribution of the temperature 
2.7 Kelvin. This radiation was predicted by Gamow as a consequence of the Big Bang. 
By the early 1970' it became clear that the CMB sky is of the dipole form (it is hotter in 
one direction and cooler in the opposite direction with the temperature difference being 
a few mili Kelvin). The dipole form can be explained by the motion our galaxy with 
regard to the rest of universe. However, at some level one expects to see irregularities, or 
anisotropies of CMB, and then it is not excluded that the generalized Planck formula will 
be appropriate for the description of CMB. 
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